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ABSTRACT 

A first-order, semianalytical theory for the long-term motion of resonant satellites is presented. 
The theory is valid for all eccentricities and inclinations and for ail commensurability ratios. 
The method aUows the inclusion of all the zonal and tesseral harmonics as well as luni-solar 
perturbations and radiation pressure. 

The method is applied to a ^nchronous satellite including only the Jjand harmonics. 
Global, long-term solutions for this problem, eccentricity, argument of perigee, and inclina- 
tion are obtained. 

INTRODUCTION 

The method of solution presented here is based on a modified Von Zeipel method applied 
to resonant satellite systems with three degrees of freedom. The method allows any in- 
clination and eccentricity and is a first-order semianalytical theory, yielding only the long- 
term perturbations. 

The procedure first eliminates the short periodic terms numerically by a classical Von Zeipel 
averaging process. Since the averaging is performed numerically, developments in the 
eccentricity and inclination are not necessary. For the same reason, it is straightforward 
to add luni-solar perturbations and radiation pressure. 

After the first Von Zeipel transformation, we have an averaged Hamiltonian that is a tabulated 
function of the variables: argument of perigee, longitude of the node, eccentricity, and in- 
clination. The averaged Hamiltonian is obtained as a numerical table, not m analytical 
function. 

The next procedure is then to eliminate the angular variable corresponding to the resonance 
argument, or the critical argument, t . a modified Von Zeipel tiansformatiort. The principle 
here is that the averaged Hamiltonian at resonance is a minimum with respect to the critical 
argument. This principle was introduced by Hoii in 1960 for his satellite critical inclination 
theory. It was used in a slightly different form by Musen and Bailie in their satellite tesseral 
resonance theory in 1 962. During the transformation, the new Hamiltonian is set equal to 
the value of the old Hamiltonian at its minimum. 
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After these two canonical transformations, a one-degree-of-freedom system remains, and 
in that system the Hamiltonian is a constant. We are now able to analyze the new one- 
degree-of-freedom Hamiltonian and obtain all long-period and secular perturbations, without 
having to develop a %ries in eccentricity or inclination. 

The procedure and algorithm of the modifled Von Zeipel method that we are using was 
developed by G. Giacaglia in 1965 (and described in detail in 1969). He applied the method 
to resonant asteroids in a circular restricted three-body problem in 1968. It was then applied 
to the elliptic restricted problems in 1970 by Giacaglia and Nacozy. The results of the 
elliptic restricted problem are in agreement with the solutions of asteroidal motion performed 
by Schubart (1968) using methods of numerical averaging. 

The method was applied to the Huto-Neptune system by Nacozy and Diehl (1974) and the 
results agree very closely with the numerical integration performed by Williams and Benson 
(1971), who used methods of numerical averaging over a four and one-half million year 
period. 

Recently, we have adapted the method to satellite systems in resonance with the tesseral 
harmonics and have applied it to the -I] *'^22 Problem. The adaptation and application are 
presented below. 

DESCRIPTION OF THE METHOD 

The method uses the Delaunay equations of motion. The angular variables are modified 
slightly so that one variable is the mean anomaly, V^. A second variable is the critical 
argument tor the synchronous satellite, Y, , and a thud i: the argument of perigee, Yj. 

Letting the X variables refer to the conjugate momentums, and where F is the Hamiltonian, 
we have 


Xo = L-H Y„ = K 

X, = H Yj = 6 + w + n+e 

Xj=G-H Yj=w 

F = 1 /2 • L* + H + geopotential (+luni-solar forces) 

The first canonical transformation is the classical Von Zeipel transformation where the new 
Hamiltonian, F*, is defined as the average of the old Hamiltonian over the mean anomaly 
from 0 to 2jt: 


Jir 

F*(Xo*.X,*,X,*;-,Y,*.Y/) - f F(X,, X,.X,; Y^, Y,, Y,)dy^,. 
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After the averaging, we replace the X and Y variables by X* and Y* and then a new Hamil- 
tonian is obtained, free from short periodic terms. This quadrature is performed numerically, 
not analytically. The quadrature is performed for various values of X^ *, Xj •, Y^ *, and Y^ *, 
and the Hamiltonian is obtained as a tabulated function of these variables. 

The new Hamiltonian, F*, has a minimum at Y^ * near 90°. This minimum corresponds to 
the libration center for Y, (also 270°). The exact value of Yj ♦ depends on the value of Y^ ♦ 
(the perigee). The partial derivative of the Hamiltonian with respect to Yj *, at Yj * near 90°, 
is zero. 

The next canonical transformation defines a new Hamiltonian, F**, to be the value of F* at 
its minimum value: 

F** (X„**, X, *•, Xj**; -, -, Yj**) 

= F* 

Y,*3s90°(or270°). 

This canonical transformation is that used in 1960 by Hori for his critical inclination solution, 
and we are using the same principle here. Since F** does not depend on Y, **, we have: 


X,- = 


3F** 

=0. 

3Y •• 


This gives a quasi-integral for Xj •* : 

Xj •* = H** = constant = >/l~c**2 cosi"*. 


This relation between e** and i** is the component of the angular momentum in the Z 
direction. Since F** is a constant of the motion (conservative system), we now have two 
relations between three variables-eccentricity, inclination, and perigee— and we can rewrite 
the two relations as: 

F**(e,i,<»)) = constant 
Ii**(e, i) = constant 

These two functions deflne curves in an e versus w plane with H** as a parameter or curves 
in an i versus cj plane, also with H** as a parameter. This then gives us the long-term solution 
for e, i, and to. 


DETAILS OF AN APPLICATION 

The first part of the numerical results for the Jj - Jjj problem are shown in figure 1 for F* 
versus the critical argument, Yj *. Only one plot is ^own here for Yj * (or cj*) 90°, but many 
other plots were obtained for different values of Yj * (<o*). We have values of eccentricity 
ranging from 0.0 to 0.9, but only six values are shown in figure 1 . Also, only one value for 
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the parameter H** is given, corresponding to H** to 0.3. It should be noted that there is 
a minimum of F* at Y, * = 90°. This corresponds to the libration center for the synchronous 
satellite at Yj • = 90°. 

We define the new Hamiltonian to have the minimum value of the old Hamiltonian, and then 
we plot the new HamUtonian. This is shown in figure 2, for F** versus Y, **, where Yj •* goes 
from 0“ to 180°. The problem is periodic inir, so this will repeat from 1 80° to 360°. 

The plot of figure 2 may be considered as a three-dimensional plot by visualizing an eccentric- 
ity axis as perpendicular to and coming out of the plane of the paper. This will then produce 
a two-dimensional surface and it can be seen that there will be a valley in the surface centered 
atYj** = 90° and e**= 0.44. 

Since F** is a constant for a given trajectOTy, we can construct planes parallel to the e versus 
cj plane (in the three-dimensional plot), one plane for each value of the constant Hamiltonian. 
In other words, we take several contour levels and each of the levels corresponds to a constant 
Hamiltonian and hence to a certain trajectory. The contour curves give the long-term global 
solutions for eccentricity versus perigee and are shown in figure 3. Cbntour levels near the 
bottom of the valley produce libration. At the bottom of the valley is the libration center, 
and at the top of the valley there is circulation. 

Consider an ec jntricity of 0.1. We see in figure 3 that, for this eccentricity, there is cir- 
culation of the perigee by following the curve, beginning at e = 0. 1 and Y^ ** = 0. In other 
words, the e** versus o)** on that curve corresponds to an actual trajectory, the long-term 
solution to the problem. 

If we have an eccentricity of about 0.4, we have libration of the nerigee about the value of 
perigee equal to 90°. At nearly e = 0.44, there is a stationary solution, and this corresponds 
to a periodic solution. Apparently, these periodic solutions have not been found prior to 
this work for the Jj - Jjj problem. 

Using the quasi-integral relation, we can obtain corresponding plots for the inclination versus 
the perigee, and this is shown in figure 4. The libration center that was shown in figure 3 
corresponding to c = 0.44 shows up at about i** = V'0° and <o** = 90°. 

The stationary solution (libration center) that we have found and presented here has an 
eccentricity of 0.44, an inclination of 70°, and a perigee of 90°. We have found that as H** 
is varied, these stationary solutions form a family of periodic solutions. 

CONCLUSION 

Our future studies will be to trace the family of periodic (stationary) solutions. We will also 
add more zonals and tesseral harmonics, the effects of the sun and the moon, and radiation 
pressure to determine what effect they will have on the family of stationary, librating, and 
circulating solutions. 
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DISCUSSION 

VOICE: What do you mean by stable stationary solutions? 

NACOZY: I mean stability in the sense that if you place a particle on any one of the contour 
curves (corresponding to a trajectory), any slight deviation from the curve will cause the 
particle merely to move to another adjacent curve. 

VOICE: What do you mean by a <«Masi-integral? 

NACOZY: The integrals that I h ; presented are accurate only to first order and hence 
are often re"'rred to as quasi-integrals. 
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